Density functional approach for inhomogeneous star polymers 
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We propose microscopic density functional theory for inhomogeneous star polymers. Our approach 
is based on fundamental measure theory for hard spheres, and on Wertheim's first- and second-order 
perturbation theory for the interparticle connectivity. For simplicity we consider a model in which 
all the arms are of the same length, but our approach can be easily extended to the case of stars 
with arms of arbitrary lengths. 
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It has been demonstrated that density functional the- 
ory (DFT) is a versatile and powerful tool to represent 
the structural and thermodynamic properties of poly- 
meric fluids 0, SUS IS W Taking into account the 
level of physical model, DFT's of polymers can be di- 
vided into two main categories. The first category 0], Q 
involves the so-called coarse-graining procedure [7J, in 
which the degrees of freedom of monomers building the 
polymer coils are integrated out. The resulting effective 
pairwise potential between the centers of masses of two 
molecules is then used in further investigations @- An 
advantage of these models emerges from the possibility 
of application of theories of simple fluids to describe the 
polymers. 

It is obvious that coarse-grained models lose some in- 
formation, e.g. a possibility of evaluation of correlation 
functions between the monomers. From this poin t of 
view, models of the second category [jj 0, IE D3 HE D3 > 
which explicitly treat the microscopic structure of poly- 
mers seem to be superior. Several microscopic DFT ap- 
proaches for inhomogeneous chain polymers have been 
roposed in the literature. Woodward and Yethiraj 
0, developed a theory that combines weighted den- 
sity approximation, known from theories of simple fluids, 
with single-chain Monte Carlo simulations. An alterna- 
tive DFT of inhomogeneous polymer solutions was for- 
mulated by Forsman, Woodward and Freasier|ll|. Their 
theory is based on the free energy functional resulting 
from generalized Flory equation of state and was ex- 
tended to the case of inhomogeneous solutions of star 
polymers |12|. However, a very convenient from numeri- 
cal point of view approach was developed by Yu and Wu 
0. This approach is based on Rosenfelds' fundamental 
measure theory (FMT) ^| and on Wertheim's first-order 
thermodynamic perturbation theory (TPT1). The the- 
ory of Yu and Wu allows for performing quite complex 
studies because it does not require single-chain Monte 
Carlo simulations and yields a fully analytical equation 
of state. This approach |(j has been successfully applied 



to investigate adsorption, surface phase transitions and 
capillary condensation in systems involving chain parti- 
cles 0[ Ell- It was also extended to the case of 
inhomogeneous semiflexible and cyclic polyatomic fluids 
0, as well as to binary hard-rod-polymer mixtures 17]. 

A few years ago Bias and Vega proposed an ex- 
tension of the associating fluid theory for bulk systems 
involving branched chain molecules. According to their 
model, branched molecules are built of chain segments 
(arms) of tangentially bonded hard spheres connected via 
articulation vertices, each of them formed by / arms. The 
excess Hclmholtz free energy due to the chain connectiv- 
ity is separated into two contributions, one accounting 
for the formation of the articulation vertex, and a second 
one due to the formation of the arms. The first term 
has been described by the second-order thermodynamic 
perturbation theory (TPT2), whereas the formation of 
chain arms - via TPT1. The principal aim of this work 
is to generalize the bulk theory of Bias and Vega to the 
case of inhomogeneous systems. We consider the sim- 
plest case of molecules with one articulation vertex. The 
generalization is carried out by utilizing the formalism of 
Yu and Wu |f|, derived for chain polymers. 

We consider an inhomogeneous fluid composed of star 
molecules, i.e. each molecule is built of a spherical "head" 
(articulation vertex), and / arms tangentially attached to 
it. Each arm is just a chain of Mf tangentially jointed 
segments. Although the numbers Mf can be different, in 
this note we study the case in which all the arms are of 
the same length, M = Mf, so that the total number of 
segments within a molecule is N — f ■ M + 1. All the 
segments are hard spheres of diameter a. The bonding 
potential H(R) is defined so that g(R) = exp[— /3VJ,(R)] 
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(1) 

where R = (r ,{r^}) with i = 1,2, •••,/ and j — 
1,2,- ■ ■ ,M denotes the set of segment positions. The 
articulation vertex is labelled by the subscript 0. All re- 
maining segments are labelled by the superscript (spec- 
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ifying arm) and the subscript (specifying the position 
within a given arm). The grand potential of the sys- 
tem 57 is as a functional of the local density of polymers, 

n[p(R)\ = F id [p(R)}+F ex [p(R)}+ [dRp(R)(V ext (R)-n), 



r (2) 

where p is the configurational chemical potential, V ex t is 
the external potential, (3F ld [p{R)] = fdRp(R)V b (R) + 
JdR/7(R)[ln( j o(R)) — 1] is the ideal part of the free energy 
and F ex is the excess free energy. The external poten- 
tial is a sum of the potentials acting on each segment, 

V ext (R) = vo(r ) + TLi ££i )■ We further as- 

sume that the excess free energy is a functional of the 
average segment local density defined as 

/ M 



p s (r) = /Oo (r) + ^^pf(r)= / dR«5(r - r )p(R 

f M 

+ EE /dIW(r-rW)p(R), 

„' 1 A 1 J 



(3) 



i=l j=l ' 



where py' (r) and po (r) are local densities of segment "j 
within the arm i" and of the articulation segment, re- 
spectively. 

Following Yu and Wu we decompose the excess free 
energy as 

(3F ex \p s {r)} = f dr[$ HS ({n a (r)}) + d> c ({n a (r)})} , 



(4) 

where §hs results from the hard-sphere repulsion be- 
tween segments, and $c is the contribution due to the 
connectivity. Each of these contributions is a function of 
four scalar and two vector weighted densities |a, Il3| . For 
the hard sphere contribution <E>#s we use the White-Bear 
theory, see Refs. 0, |2(j for the explicit formula. 

Wertheim's perturbation theory for a bulk fluid p2] | 
can be naturally incorporated into the DFT framework 
|2l| . The generalization for inhomogeneous star polymer 
systems is represented by the expression 0, |23] 

<£ c ({n a (r)}) = $ arm ({n a (r)}) + $ art ({n tt (r)}) , (5) 

where $ arm and Q art are the contributions due to the 
formation of chains within consecutive arms and due to 
the formation of the articulation vertex. The equation 
for <j> arm follows from the theory of Yu and Wu [6j 



$ arm ({n a (r)}) 



1 + f-N 
N 



n (\n[y H s(<J, {n a (r)})}, 



(6) 

where £ = 1 — ny2 • ny^A?^) 2 , and the contact value of 
the hard-sphere cavity function, dhs{v)i results from the 
Carnahan-Starling equation of state, cf. Eq.(18) of Rcf. 
0. Free energy density <& art is obtained by generalizing 
the theory of Bias and Vega [Tll^ 

® art ({n a (r)}) 



® a T%i({n a (v)}) + ^ T2 ({n a (T)}), 
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FIG. 1: The average segment density profiles of star poly- 
mers (/ = 3, M = 5) evaluated for the bulk segment packing 
fraction r/ a b = 0.1,0.2 and 0.3. Symbols represent computer 
simulations |24|, dashed lines denote DFT results obtained 
using TPT1 and dots - DFT results evaluated using TPT1 
and TPT2 contributions. 



where $5? T1 and ^pj^ represent the first 

$^ T1 ({n Q (r)}) = -ln CHy HS (a,{n a (r)})] , (8) 
and the second-order perturbation terms jTsL l22T | 



cb art 



In Vl + 4A (9) 

(l + VTTW+ 1 - (i - VTT4A) /+1 



In- 



2/H 



(7) 



In the above A depends on the number of arms and its 
evaluation requires the knowledge of the /-body corre- 
lation function for / spheres in contact. In the case of 
/ = 3 the application of the superposition approxima- 
tion yields A = (1 + a nz + 6n|)/(l — n^) 3 — 1, where a 
and b are constant that depend on the angles between 
the arms attached to the articulation vertex 0, |2^, |23| . 
In the case of bulk fluids ri3 is just the packing fraction. 
Approximation proposed for inhomogeneous system re- 
lies on substitution of the bulk packing fraction by the 
weighted density 77,3. Note that this approximation is 
not unique. One can follow the ideas of Yu and Wupof 
and propose an approximation involving, besides scalar, 
also vector weighted densities. In this work, however, we 
decided to apply as simple expression, as possible. 

Within the TPT1 approach the bulk thermodynamic 
properties of the star polymers are the same as the prop- 
erties of chains built of the same number of segments. 
This is because the first-order perturbation free energy 
takes into account only the number of segment connec- 
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tions and neglects polymer's topology. The latter is in- 
cluded within the TPT2 approach, cf. Eq.(15) of Ref. 
|18| . However, the identity of the bulk thermodynamic 
properties within the TPT1 theory does not imply that 
the structure of inhomogeneous fluids of star polymers 
and of chains with the same number of segments, that 
results from DFT, is identical. 
Minimization of J2J yields 



z/o 



P( R ) = 5(R) ex P 



/ M 



/?/i-/3Ao(r )-/3^X>i (r 



1=1 j=l 



(10) 

where Af (rf) = 5F ex /5p s {vf )+vf{rf) and A (r ) = 
SF ex /5p s (ro) + va(ro). For systems with the density dis- 
tribution varying only in the z direction Eqs. and l|l(J|) 
lead to the following expressions for the segment density 
profiles 

Po(z ) = exp( / 3 / i) 7 o(zo) [G M+1 (z )] f (11) 



and 



pf{zf ) ^^{^'{z^G^-^z^G^^') , 

(12) 

where 7 j l) (z) = exp[-pxf ] (z)]; l0 (z) = ^{z). The 
functions G l (z) are defined by recurrence relation 6] 



(P(z)= I d*' 7 jV) 



(a - \z - z'\) , 



2a 



G 3 {z ) , (13) 



for j = 2, • • • , M with G i {z) = 1. In the above 9 is the 
unit-step function. The functions &{z), however, are 
given by 



•■■•2, . . _ --'/' K7 l z z l ' ) -[G' M+1 (z')] /_1 



G 2 (z) = / dz' 7o (*') : 



2(7 



for j ' = 2 and 

G J (z) = / d*' 7 f ( 



(14) 



'-w^fcdizf^-V), (is) 

2ct 



for j > 2. The equations given above are valid for the 
stars with arms of identical length. In such a case the 
profiles Pj (z) are independent of the arm index i. A 
generalization of the theory to the case of stars with 
arms of different length is straightforward. For exam- 
ple, the integrand in the last equation would involve a 
product of the functions G Ml+1 (z')G Al2+1 {z') . . ., instead 
of [G M+1 (z')}'~ 1 (here Mi's abbreviate the number of 
segments within consecutive arms). As a simple appli- 
cation of the theory we calculate density profiles of star 
molecules built of hard-sphere segments near a hard wall, 
located at z = 0. The solutions of the density profile 
equations were obtained by using the standard iterational 
procedure. 

In Fig. 1 we compare the average segment density pro- 
files resulting from theory with computer simulations p4| 




FIG. 2: Upper panels: The average segment density pro- 
files of star polymers (dash-dotted lines) and of chains (solid 
lines) calculated for the total number of segments N = 13 
and N — 61. The results for N — 61 are shifted up by 0.2. 
Lower panels: The segment density profiles of the articulation 
segment (circles) and of its nearest neighbor (squares) of star 
polymers and of the first (solid lines) and the second (dashed 
lines) segment of chains In left panels p* b — 0.2, whereas 
p* sb — 0.6 in right panels. 



for star polymers built of / = 3 arms, each composed of 
M = 5 segments. The calculations were carried out for 
bulk segment packing fractions rj s i> — irp s t,a 3 /6 = 0.1, 0.2 
and 0.3, where p s {, is the bulk average segment density. 
The density profiles in Fig. 1 have been normalized by 
the bulk density p s b. For 77^ = 0.2 and 0.3 we show two 
sets of the DFT results. The first one has been evalu- 
ated employing the TPT1 approach (i.e. the term given 
in Eq.|H|has been neglected), whereas the second set was 
obtained using TPT1 and TPT2 contributions to the free 
energy. The differences between the local densities result- 
ing from these two approximations are small and occur 
only within the region adjacent to the wall. The TPT2 
contribution leads to smaller contact values of the aver- 
age segment local density. This effect is quite obvious, 
because the TPT2 correction lowers the pressure. The 
agreement between theoretical predictions and computer 
simulations is reasonable. 

Figure 2 compares the density profiles for three-armed 
star polymers (resulting from the TPT1 approach) with 
the profiles of chain polymers built of the same number 
of segments obtained from the approach of Yu and Wu 
The results are for bulk average segment densities 
P* s b = Psb " 3 = 0.2 and 0.6 and for two model systems 
with different total number of segments N = 13 and 
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FIG. 3: (a) The average segment density of star polymers 
with / = 3 and M = 20 (dashed line) and with / = 4 and 
A/ = 15 (solid line) for two bulk densities, p* b — 0.2 and p* sb = 
0.4. (b) The segment density profiles of "heads" (solid and 
dash-dotted lines) and it's first neighbor (dashed and dotted 
lines) for the same models. The bulk segment density is p sb = 
0.4. 

N = 61. For N — 13 each star polymer arm is composed 
of M = 4 segments, whereas for N = 61 each arm con- 
sists of M = 20 segments. In the upper panels we com- 
pare the average segment densities normalized to unity. 
We find that for higher bulk density p* b — 0.6 (the upper 
right-hand side panel) the local densities p s {z) of chains 
and stars are quite similar. The profiles are dominated 
by packing effects. Larger differences between the star 
and chain polymer profiles are visible at lower density, 
p* b = 0.2, cf. the upper left-hand side panel. Note that 
the contact values of p s (z) for star and chain polymers 
are identical in our TPT1 approach. Lower panels of Fig. 



2 show the density profiles of selected segments for the 
same systems. We plot here the profiles of "heads" (in 
the case of chains the profiles of the first segment) and the 
profiles of the segments attached directly to the "heads" . 
The differences between the profiles for the chain and 
star polymers are now more pronounced, especially for 
the profiles of the "heads". We have also inspected the 
profiles for the segments that are topologically more dis- 
tanced from the head and have found that the difference 
between them becomes gradually smaller. 

Finally Fig. 3 presents the profiles of the stars built of 
the same number of segments, but having different num- 
ber of arms. We have considered the models with / = 3, 
M = 20 and with / = 4, M = 15. It is not surprising 
that the difference between the average segment density 
profiles (cf. Fig. 3a) is now less pronounced than in the 
case of the profiles displayed in Fig. 2, because the dif- 
ferences in the topology of the particles are now smaller. 
However, the differences between the individual segment 
density profiles still persist, cf. Fig. 3b, where we show 
the profiles of "heads" and the segments directly attached 
to heads. 

In conclusion, in this work we propose density func- 
tional theory for inhomogeneous star polymers. Al- 
though the theory is written down for the case of arms 
composed of identical numbers of segments, its general- 
ization for stars with arms of different length is straight- 
forward. Several further extensions are also possible. In 
particular it would be of interest to consider cases of 
physically different "head" and "arm" segments in or- 
der to mimic the systems involving surfactants. Some of 
these topics are already under study in our laboratory. 

This work was supported by EU under a TOK contract 
No. 509249. We are grateful to prof. O. Pizio for early 
discussions and bringing Ref. [lj| to our attention. 
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